Nonlinear states and nonlinear tunneling 
in a potential well 
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A nonlinear Schrodinger model in a square well and managed nonlinearity is shown to possess 
nonlinear states as continuous extensions of the linear levels. The solutions are remarkably stable up 
to a threshold amplitude where a soliton is emitted and propagates outside the well. The analytic 
expression of the threshold is given in terms of the well size for each level. This process of nonlinear- 
tunneling results from an instability of the evanescent wave inside the walls and can find experimental 
realization in a proposed nonlinear fiber Bragg gratings resonator. 



Introduction. We consider the nonlinear Schrodinger 
equation (NLS) for the wave function ip(x, t) trapped in 
a square well, in reduced units, according to the following 
model 



|z| > L 
\x\ < L 



itpt + ipxx + M V = Vtp, 
iipt + ipxx - MV = 0. 



(1) 



The positive constant V is the well height and 2L its 
width. The nonlinearity is managed such as to be repul- 
sive inside the well and attractive outside. 

The nonlinear states for NLS with a definite sign of 
the nonlinearity has been derived in Q with interesting 
solutions having no counterparts in the linear limit. Us- 
ing nonlinearity management as in Q actually allows to 
obtain a system whose nonlinear states solutions possess 
remarkable properties: i) they uniformly tend to the lin- 
ear eigenstates in the small amplitude limit, ii) they are 
stable solutions for amplitudes below an explicit thresh- 
old, iii) above the threshold an instability generates gap 
solitons, emitted outside the well, hence realizing a clas- 
sical nonlinear tunneling process. 




FIG. 1: Example of a physical situation where the model (0 
applies where t is the propagation direction and V is propor- 
tionnal to the dielectric constant difference eo — ei. 

This problem has promising application to nonlinear 
optics where it applies to the propagation in the direc- 
tion t, and tranverse modulation along x, of the enve- 
lope of a laser beam in a fiber Bragg grating (in Kerr 
regime) 2] when different types of fibers are used: self- 
focusing nonlinearity with dielectric constant e\ outside 
[— L,+L], self-defocusing with eo > e± inside. The peri- 
odic medium extending in \x\ > L acts as a Bragg mir- 
ror for the tranverse modulation and the resulting Bragg 



grating resonator is sketched on figure ^ The arrows 
show the injected radiation. 

The domain of spatial solitons in fiber Bragg gratings 
is rich of recent results, see e.g. 0,with n i ce exper- 
imental evidence of soliton formation |j| and recent ex- 
periments which demonstrated the existence of discrete 
stationnary gap solitons @- On the theoretical side, an 
interesting proposal of Q is to generate discrete gap soli- 
tons by boundary driving a Bragg grating above the cut- 
off. The aproach uses the theory of nonlinear supratrans- 
mission |8j extended to the case when the forbidden band 
results from the discrete nature of the medium. 

In the quasi-continous case (large number of fibers and 
weak coupling), the device of fig Q] will be shown to be a 
means to generate gap solitons by continuous wave (cw) 
input radiation at a flux intensity given explicitely in 
terms of the dimensions of the potential well depth V (or 
the dielectric constant varaition eo — ei). As an illustra- 
tion, figEl shows the evolution of the fundamental mode 
at the threshold for which the input envelope is indeed 
cw-like and experiences an instability generating the gap 
soliton. The threshold vaues are explicitely calculated for 




FIG. 2: Profile of F — \ip\ for the fundamental 0-mode initial 
datum at threshold amplitude in the case V = 1 and L — 4. 

the entire set of nonlinear states. Finally we demonstrate 
that the linear small amplitude limit maps continuously 
the nonlinear states to the linear ones. 

We shall need to refer to the linear eigenstates that 
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read in the odd case 

\x\ < L : 
\x\ > L : 
k 2 = Vsm(kL) , 

and in the even case 

\x\ < L : 
\x\>L : 
k 2 = Vcos(kL) , 



i/j = Ae- iut sin(kx) , 

i/j = Ae- iut sm{kL)e~ K ^ 

tan(fcL) < . 



(2) 



ip = A, e~ lult cos(kx) , 

^ = Ae- iut cos{kL)e- K ^- L ^ , 

tan(fcL) > , (3) 



for Co = k 2 < V and k 2 — V — lu, and for arbitrary 
amplitude A. 

Nonlinear states. As learned from y|, the main tool 
to derive the nonlinear states is to connect a periodic 
solution inside the well to a static one-soliton tail out- 
side. We make use of the fundamental solutions of NLS 
in terms of Jacobi elliptic functions for both symmetric 
(even) and antisymmetric (odd) cases @. 

First, outside the potential well we have the common 
soliton tails given by (still with n 2 = V — ui and ui < V) 



\x\ > L, ip 



cosh[«;(|a;| — L + d)] 



(4) 



which replaces the two evanescent waves of the linear 
case. Note that the function (—tp) is also solution, a 
property which must be used to conveniently connect the 
solution inside to the outside tails. 

Inside the well we use the basic solution 



< L, V = LLAe~ iuJt sn(A'x + x ,(J,). 



(5) 



Here and in the following the amplitude A is real- valued 
and positive and we define the new constant A' by A = 
A'\[2. By convenient choices of the constant Xq we shall 
obtain both the odd and even solutions in the well. 

This function ip is a solution of if (necessary con- 
dition) 



c 2 = ±A 2 (l + M 2 ), 



(6) 



which links the frequency lo to the modulus /i of the ellip- 
tic functions. The requirement u> < V for a bound state 
implies that the modulus /z cannot exceed a maximum 
value n m 



/./, 



< (2V/A 2 - l) 1 / 2 = 



Mr, 



(7) 



We discovered that it is necessary to work with ex- 
pressions of the solution for moduli fj, > 1, performed by 
means of the identity /xsn(u, /i) — sn(/iu, l//i). The so- 
lutions for fi < 1 can be easily written down but either 
they do not contribute to the spectrum or they have their 
exact conterpart with fj, > 1. 



Odd solutions. In the range /j £ [1, /i m ] we define the 
following solution inside the well 

Id < L, = Ae~ lu}t sn(/iA'x, -) . (8) 

In order to obtain sufficient conditions for the set JSJ) (@J 
to be a solution of we need to ensure continuity of 
the solution and its derivative in x — ±L. This provides 
the admissible discrete set of moduli fi which, after some 
algebra, must solve 

^ 2 cn 2 (b,-)dn 2 (b.i) = 



sn 2 (b. 



1 



2V 
A2 



(l+/i 2 )- sn 2 (b,-) 



where the function b is defined as 



HA'L. 



(9) 



(10) 



Then the shift d of the tail position for each solution /i 
of the above equation, is given by 



A sn(b,-) = — —, —r . 
fj, cosh(Kd) 



(11) 



As for the linear case we must add a consistency con- 
dition which ensures the continuity of the derivatives. It 
is just a matter of careful reading of all possible cases to 
obtain the condition 



sn(b, -) cn(b, -) dn(b, -) < 0, 
A* A 4 A* 

for the solution fi of © to be acceptable. 

Even solutions. Still with \x e [l,/i m ] we define 



(12) 



\x\ < L, tp = At 



~ iut sn( At A'x + K(-),-) 
A* A* 



(13) 



obtained from JSJ) by the translation xo = K(1//j,)//j. 
where K is the complete elliptic integral of the first kind. 
In that case the continuity conditions give the admissible 
set of moduli as solutions of 10 with now 



b = fiA'L + K(-) 
A* 



(14) 



and the consequent definition l|llfl of d. The consistency 
condition H12fl still holds here with the above function b. 

Note that in both expressions JSJ and , the param- 
eter A does represent the maximum value of the ampli- 
tude of the solution, reached at x — for the even states 
(0 and at x = K(\j\i) for the odd states ©• 

Example. The figure shows the solutions of the 
equations © for V = 1 and L = 4 in the odd and even 
cases, for which the associated linear problem possess 3 
eigenstates (dashed lines) . The 3 curves show the depen- 
dency of the eigenvalue of each nonlinear extension of the 
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FIG. 3: Dependence of the eigenvalues ui in terms of the am- 
plitude for V = 1 and L — 4. Crosses indicate the values of 
amplitude threshold and the dashed curve is the prediction of 
global threshold (11911 . 



linear levels in terms of the amplitude A. These eigen- 
values are given by the expression © for each solution 
/i of (JSJ, in both cases l|10(l and l|14f) . equations that are 
solved numerically (like in the linear case). 

The next figure 0] displays the three states for V = 1 in 
the linear and nonlinear case corresponding to particular 
choices of amplitudes and frequencies as indicated on the 
graphs. These are the plots of the solutions © and (jT3|l . 
compared to the linear eigenfunctions © and J2J|. 



LINEAR 



NONLINEAR 




-0.5 - A=0.5 



FIG. 4: Examples of linear and nonlinear states, for V = 
1, with amplitudes and frequencies as indicated. The same 
vertical scale is used for linear and nonlinear figures. 

There are two fundamental properties of this "nonlin- 
ear spectrum" that we demonstrate hereafter. It is first 
the evident natural property to reach the linear spectrum 
for vanishing amplitudes. This explains in particular why 



we have the same number of states as the linear levels. 
Second all three curves are seen to stop at some threshold 
amplitude which can be exactly computed, as shown by 
the dashed curve on figure |21 and the crosses which are 
the predicted thresholds. 

Linear limit. To demonstrate in general that the non- 
linear spectrum goes to the linear one in the limit A — ► 0, 
it is crucial to evaluate the limits keeping the product Afj, 
finite. Actually defining 

k = -j=n, (15) 

we can easily obtain (A = A'y2) 

Asn(fiA'x, 1/n) ~ q Asin(kx), (16) 

which proves that the nonlinear odd states JHJ tend to 
the linear ones in the limit A — > 0. The same property 
holds naturally for the nonlinear even states (|13f) by using 
sn(u+K(i/), v) = cd(u, v) and cd(u, 0) = cos(u) together 
with l|16l) above. 

These results show in particular that the "nonlinear 
wave number" is the quantity k — A' fi and thus that the 
relation |JBJ must be read oj 2 — k 2 + A 2 /2 that tends 
to the linear dispersion relation u> 2 — k 2 in the small 
amplitude limit. 

Nonlinear tunneling. More interesting is the existence 
of a threshold amplitude (or population threshold) be- 
yond which gap solitons are emitted outside the well 
hence realizing a classical nonlinear tunneling. This is 
the result of an instability, as described in [l(j, which 
takes place as soon as the amplitude A is such that the 
soliton tails Jlj reaches its maximum value in |x| = L, 
i.e. d = 0. 

In order to derive the analytic expression of the thresh- 
olds postions {A s ,u) s } for each branch lu(A), we express 
the continuity conditions in \x\ = L in the case d = 0. 
The threshold is thus obtained as a particular solution of 
the equation for the state © for which both sides van- 
ish altogether (successively for the odd and even states), 
namely 



cn 2 (b,-)dn 2 (b,i) 



0, 



--C 

A 2 1 



sn 2 (b, -) 



0. 



(17) 
(18) 



where b is given by IjlOjl for the odd solutions and by l|14|l 
for the even ones. 

It is then a matter of algebraic manipulations to 
demonstrate that the solutions {A s , fi s } of the above 
equations globally satisfy 



At 



2V 



2(V-w.), 



where /i s is obtained by solving for the odd case 



sn 2 (L/i<t 



V 



1 



M 2 + 2' n' 



1, 



(19) 



(20) 
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and for the even case 

sn 2 (W^— +K(i),i) = l. (21) 

y n 2 + 2 n fi 

Note that the solutions {A s , fi s } of these equation must 
obey requirement (JZJ) which is checked a posteriori. 

The dashed curve of fig. |3|is the plot of JTJJ for V = 1 
and L = 4, and the crosses are obtained by solving l|20|) 
and H21fl numerically 

The above pocedure does not furnishes the threshold 
corresponding to the fundamental level. Indeed it misses 
the particular solution 

Ms = l, A l = lv, (22) 

which is effectively a solution of (|17(l(|18|l from the prop- 
erty 

lim(sn(a + K(-),-)) = l, (23) 

for any real-valued a. In that case the "dispersion re- 
lation" © provides uo s — A 2 . Note that the threshold 
amplitude A s does not depend on the width 2L of the 
well. 

This solution is particularly interesting in view of ap- 
plications to the fiber guide grating of figH Indeed in 
that case the corresponding expression l(T3*|l is the con- 
stant amplitude field ip — A s e~ lbJst (which is trivialy a 
solution of (3). Injecting then in the medium (as indi- 
cated by the arrows on fig^l a cw-laser beam of (nor- 
malized) intensity flux A 2 constant along the tranverse 
x-direction for x € [— L,+L], one would generate a gap 
soliton propagating outside the well, as displayed in fig. 
[5] Such simulation can be reproduced ad libidum for 
input data constant in [— L, +L] and exponentially van- 
ishing outside. As soon as the amplitude exceeds the 
threshold A s , one obtains tunelling by soliton emission. 

It is clear that the present formalism has been devel- 
opped for continuous envelopes t/j(x,t), which requires a 
large number of fiber guides and weak transverse cou- 
pling. An interesting extension is then the discrete situ- 
ation with the open question of the discrete states solu- 
tions. 



Perspectives and conclusion. We have so far demon- 
strated the existence of "nonlinear eigenstates", exact 
solutions of the NLS model JIJ, that are remarkably sta- 
ble (numerical simulations did not show any deviation 
from the exact expressions up to times when numerical 
errors become sensible, i.e. 10 3 to 10 4 for our scheme) 
as soon as their amplitudes do not exceed a threshold 
explicitely evaluated in terms of the well height for each 
mode. These states reproduce exactly, in the small ampli- 
tude limit, the usual eigenstates of the Shrodinger equa- 
tion in a potential well. 

Above the threshold amplitude, a "nonlinear tunnel- 
ing" (macroscopic, classical) occurs by the emission of 
gap solitons outside the well. This generic property sug- 
gests a "gap soliton generator" by the device sketched in 
figd more especially as the fundamental level, close to 
the threshold, behaves as a cw-field. 

The approach applies straightforwardly to a model 
where the inside of the well would obey the linear 
Schrodinger equation. This in view of applications to 
a resonator device with Bragg medium playing the role 
of the mirrors, a subject under study now. 

Another natural question is the properties of the non- 
linear states when the inside of the well presents an at- 
tractive (focusing) nonlincarity (as the outside medium) . 
In that case, similar states are defined in terms of Ja- 
cobi cn-functions but we have observed that they 
experience modulational instability before reaching the 
threshold amplitude. In the other case when the whole 
medium is defocusing, the outside tails are of cosech-type 
and thus can match any amplitude of the nonlinear state, 
hence there is no threshold and no gap soliton generation. 
Thus the management of the nonlincarity chosen in Q is 
a means for sabilizing the nonlinear states and ensuring 
a threshold for soliton generation. 

We report to future studies the important question of 
the analytical proof of the stability of those states, includ- 
ing the mathematical derivation of the instability crite- 
rion above the threshold amplitude. Note finally that, 
in the threshold regions, the states are degenerated: two 
solutions coexist at the same frequency with slightly dif- 
ferent amplitudes as seen on fig EH 
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